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Abstract 

Splines can be constructed by convolving the indicator function of 
a cell whose shifts tessellate K". This paper presents simple, non- 
algebraic criteria that imply that, for regular shift-invariant tessella- 
tions, only a small subset of such spline families yield nested spaces: 
primarily the well-known tensor-product and box splines. Among the 
many non-refinable constructions are hex-splines and their generaliza- 
tion to the Voronoi cells of non-Cartesian root lattices. 

1 Introduction 

Univariate uniform B-splines can be defined by repeated convolution, start- 
ing with the indicator functional of the intervals or cells delineated by knots. 
This construction implies local support and delivers a number of desirable 
properties (see [dB78, dB87]) that have made B-splines the representation of 
choice in modeling and analysis. In particular, B-splines are reflnable. That 
is, they can be exactly represented as a linear combinations of B-splines with 
a finer knot sequence. Refinability is a key ingredient of multi-resolution and 
adaptive and sparse representation of data. Refinability also guarantees 
monotone decay of error when shrinking the intervals. 

By tensoring univariate B-splines, we can obtain splines on Cartesian 
grids in any dimension. Box-splines [dHR93j generalize tensoring by allow- 
ing convolution in directions other than orthogonal ones. As a prominent 
example in two variables, the linear 3-direction box-spline consists of linear 
pieces over each of six equilateral triangles surrounding one vertex. Shifts of 
this 'hat function' on an equilateral triangulation sum to one. Convolution 
of the hat function with itself results in a twice continuously differentiable 
function of degree 4; and m-fold convolution yields a function of degree 
3m — 2 with smoothness C 2m . Since this progression skips odd orders of 
smoothness, van der Ville et al. |vBU + 04] proposed to directly convolve 
the indicator function of the hexagon and build splines customized to the 



: An indicator function takes on the value one on the interval but is zero otherwise. 
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(a) natural and man-made hex-tilings 



(b) half-scaled hex-tiling 



Figure 1: Hexagonal tessellations. (a) basalt formation and tiles 
(http:/ /en.wikipedia.org/wiki/{Basalt,Hexagonal_tiling}) (b) Non-nesting 
of the hex partition in Example 1. 

hexagonal tessellation of the plane (cf . Fig. [Th.) . They went on to show that 
the resulting hex-splines share a number of desirable properties familiar from 
box-splines. But the authors did not settle whether the splines were refinable 
[vdVUlO] . i.e. whether hex-splines of the given hexagonal tessellation T can 
be represented as linear combinations of hex-splines based on a scaled-down 
hexagonal tessellation, say ^T. Generalizing the analysis of hex-splines, 

• this paper presents simple non-algebraic criteria necessary for regular 
shift-invariant tessellations to admit refinable indicator functions. 

For example, such a tessellation must contain, for every cell facet /, the 
plane through /. Therefore, requiring refinability, even of just the constant 
spline, strongly restricts allowable tessellations. 

• In contrast to tensor-product and box splines, we show that hex-splines 
and similar constructions can only be scaled, but not refined: scaled 
hex-spline spaces are not nested. 

• The analysis extends to overcomplete families (superpositions) of spline 
spaces. 

The following example illustrates how non-refinability leads to loss of mono- 
tonicity of the approximation error under scaling: for one or more steps 
halving the scale can increase the error. By contrast, nested spaces guaran- 
tee monotonically decreasing error. 

Example 1 Let T~L l be the space of indicator functions over a regular tes- 
sellation by hexagons of diameter 2~ l and such that, at each level of scaling, 
the origin is the center of one hexagon. Denote by H the indicator function 
in %° whose support hexagon is centered at the origin. H 1 does not contain 
a linear combination of functions that can replicate H since the supports 
of the six relevant scaled indicator functions are bisected by the boundary 
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of the support of H (see Fig. [T^). Correspondingly, the L 2 approximation 
error to H from H 1 is ^A 1 > where A 1 is the area of the hexagon with 
diameter |. Since the error from T-L° is by construction zero, the scaling by 
1/2 has increased the error. By carefully adding to H an increasing number 
of scaled-down copies, small increases in the L 2 error can be distributed over 
multiple consecutive steps. □ 

Overview. Section [2] reviews tessellations induced by lattices, hex- 
splines and their generalizations. Section exhibits two non-algebraic crite- 
ria, chosen for their simplicity, for testing whether a tessellation can support 
a refinable space of splines that are constructed by convolution of indica- 
tor functions of its cells. Section H] extends this investigation to a multiple 
covering of M. n by distinct families of indicator functions. 



2 Splines from lattice Voronoi cells 

A n-dimensional lattice is a discrete subgroup of full rank in a n-dimensional 
Euclidean vector space. Alternatively, such a lattice may be viewed as induc- 
ing a tessellation of space into identical cells without n-dimensional overlap^. 
The tessellation is then generated by the translational shifts of one cell. For 
example, lattice points can serve as sites of Voronoi cells. The Euclidean 
plane admits three highly symmetric shift-invariant tessellations: partition 
into equilateral triangles, squares, or hexagons respectively. Repeated con- 
volution starting with the indicator function of any of these polygons yields 
spline functions of local support and increasing degree. The regular par- 
tition into squares gives rise to uniform tensor-product B-splines and the 
regular triangulation and its hexagonal dual to box splines. 

An interesting additional type of spline arises from convolving the indi- 
cator function H of the hexagon with itself. Such hex-splines, a family of 
C fc_1 splines supported on a local k + 1-neighborhood, were developed and 
analyzed by van De Ville et al. |vBU+04] . That paper compares hex-splines 



to tensor-product splines and uses the Fourier transform of hex-splines to 
derive, for low frequencies, the 1? approximation order, as a combination 
of the projection into the hex-spline space and a quasi-interpolation er- 
ror. [CvB05j derived quasi-interpolation formulas and showed promising 
results when applying hex-splines to the reconstruction of images (see also 
|CvU06l ICvOTl EvQ8]). Van De Ville et al. [vBU+04] . also observed that 



hexagons are Voronoi cells of a lattice and that the cell can be split into three 
quadrilaterals, using one of two choices of the central split. Thus H can be 
split into three constant box splines whose mixed convolution yields higher- 
order splines [Kim08b, ME10]. However, while box-splines are refinable, we 
will see that hex-splines are not refinable in a shift-invariant way. 



2 A common convention is to define the cells to be half-open sets so that they do not 
overlap on facets, but nevertheless cover. 
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3 Refinability constraints 



We consider a polyhedral tessellation T of W 1 into unpartitioned re-dimensional 
units, called cells, that are bounded by a finite number of n — 1-dimensional 
facets. We denote by x(T) the space of indicator functions of the cells of T 
and by xC^ 1 ) the space of indicator functions on some scaled-down copy T 1 
of T. The space x(T) is refinable if each indicator function in x(T) can be 
represented as a linear combination of functions in ^(T 1 ). 

Establishing whether a tessellation T admits a refinable space of indica- 
tor functions therefore requires proving the existence of weights such that a 
linear combination of elements in x^ 1 ) with these weights reproduces each 
element in x(T). Proposition [TJ below provides a much simpler necessary 
condition that avoids such algebraic analysis. While our focus is on shift- 
invariant tessellations, Proposition [1] applies more generally and also to cell 
boundaries of co-dimension greater than 1. Its proof uses the notion of a 
c 1 G T 1 straddling a facet of a cell c G T. A cell c 1 straddles a facet / of c 
if there exists a point p on /, a unit vector n normal to / at p and e > 
such that both p + en G c 1 and p — en G c . 

Proposition 1 Lei T be a polyhedral tessellation of W 1 and T 1 zis scaled- 
down copy. Then x(T) is refinable only if every facet of T is the union of 
facets of T 1 . 

Proof Assume that a facet / of a cell c in T is not a union of facets of 
T 1 . Then, since T 1 is a tessellation, some cell c 1 of T 1 must straddle /. Let 
H 1 G xl^ 1 ) be the indicator function of c 1 and H the indicator function 
of c. Then, in order to reproduce the unit step of H across /, H 1 must 
simultaneously take on both the value and the value 1. J J 



Translation-invariant or shift-invariant tessellations are a special case 
of transitive tilings where every cell can be mapped to every other cell by 
translation, without rotation. 

Proposition 2 If T is a shift-invariant tessellation, x(T) is refinable only 
if, T contains, for each facet f , the hyperplane through f . 

Proof The coarser-scaled copies of T contain enlarged copies of every facet 
in T. By Proposition [1] these copies must be a union of facets of T. There- 
fore a shifted copy of every facet is strictly contained in the interior of and 
so extended by some coarser facet. Shift-invariance then implies that ev- 
ery facet / lies strictly inside such an extension. Ever coarser tessellations 
provide a sequence of extensions of / in any direction by any amount. || 



By inspection of the three regular tessellations of the plane, only the 
Cartesian grid and the uniform triangulation satisfy Proposition [21 but not 
the partition into hexagons. 
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Corollary 1 Hex splines are not refinable. 

We can generalize this observation by simplifying the inspection criterion. 




Figure 2: A pair of facets fi, fi of a cell c meet with an obtuse angle: 
the outward pointing normals ni, ri2 (dashed) have a strictly positive inner 
product. 



We say that two abutting facets f\ and fi of a cell c meet with an obtuse 
angle if, for i = 1,2, there exist unit vectors rij orthogonal to fi and outward 
pointing so that ni • ri2 > 0. 

Proposition 3 Let T be a tessellation ofW 1 by shifts of one polyhedral cell 
c. If two facets f\ and $2 of c meet with an obtuse angle and if d , the 
reflection of c across fi, is a cell of T then x(T) is not refinable. 

Proof Assume x(T) is refinable under the given conditions. Let d be the 
reflection of c across (the plane through) fi- Since d must not overlap c, 
obtuse angles exceeding tt, such as the reentrant corner of an L-shaped cell, 
cannot occur in c. Denote by /{ the reflection of f\ across fi and by e the 
common intersection of f\ , f% and /{ (see Fig. [2]) . 

Within d , by reflection, the facets fi and f[ meet at e with an obtuse 
angle. By Proposition [2] the extension F\ of f\ lies in T. Since the outward- 
pointing normal of fi with respect to d is — n2, fi and F\ meet at e with an 
acute angle. Therefore F\ extends f\ into and splits d . This contradicts the 
definition of a cell as an unpartitioned unit and hence the initial assumption. 

The next Proposition [3] allows us to quickly decide which of the (symmet- 
ric crystallographic) root lattices A n , A* n , B n , V n , P*, £j, j = 6, 7, 8 [CS98] 
are suitable for building refinable splines by convolution of their nearest- 
neighbor (Voronoi) cells. 

Corollary 2 Splines obtained by convolving the Voronoi cell of a non-Cartesian 
crystallographic root lattice are not refinable. 
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Proof We test whether the Voronoi cells of the root lattices contain a pair 
of abutting faces that meet with an obtuse angle. We may assume that one 
Voronoi site (cell center) is at the origin. By definition of a Voronoi cell, the 
position vectors ni and 112 of two adjacent nearest neighbors, as identified by 
their root system, are normal to the corresponding abutting bisector facets. 
Therefore these facets meet with an obtuse angle if ni • ri2 > 0. 

The An lattice is traditionally defined via an embedding in M n+1 , n > 1. 
More convenient for our purpose is the alternative geometric construction 
in M n via the n x n generator matrix A n := I n + ^J n of Theorem 1 of 
[KPlOj . Here I n is the identity matrix, J n the n x n matrix of ones and 
c n := _ Denoting the ith coordinate vector by ej, we choose ei 



and ei + e2 on the Cartesian grid, and map them via A n to the nearest An 
neighbors of the origin. The inner product of the images of ei and ei + e2 
is 

A w e x - A n ( ei + e 2 ) = + n+ " = -(n + v^TT- 1) > 0. 



n n 

1+- 



For the A* n lattice, the computation is identical except that c n : 
The inner product is w ^ n 1 +1 ^ {n 2 — 2n-2 + 2 v / ^TI)) > 0. 



For the T> n lattice, defined in n > 3 dimensions, the generator matrix is 



D r 



T 1 -p n_1 

i-n-l G n-1 

-i* -1 

Jn-l 1 



(see e.g. Section 7 of jKPllj ) and 

D n ei -D n (ei +e 2 ) = 3 > 0. 

Since D~* is the generator of T> n , the inner product for D n is 2. 

For B n , the Cartesian cube lattice has an inner product of identifying 
its uniform tensor-product B-spline constructions as potentially refinable 
(which indeed they are). On the other hand, splitting each cube by adding 
the diagonal directions of the full root system |Kim08a] yields the inner 
product ei • j = 1. 

For £q, we select the root vectors (1,1,0,0,0,0) and (1, 1, 1, 1, 1, y/S)/2 
with inner product 1. For £7, we select the root vectors (1, 1, 0, 0, 0, 0, 0) and 
(l,l,l,l,l,l,v / 2)/2 with inner product 1. For £g, we select the root vectors 
(1,1,0,0,0,0,0,0) and jg/2 with inner product 1. || 

The equilateral triangulation in R 2 is dual to the 'honeycomb lattice' 
which is not a standard lattice. The equilateral triangulation yields an 
inner product of ^ compatible with refinability and indeed plays host to 
the refinable 'hat' function. 



4 Overcomplete spaces 

Since the evaluation of hex-splines by convolving three families of box splines 
already makes use of a large number of terms, it is reasonable to investigate 
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whether superposition of several families of hex-splines are refinable as a 
family. That is, we consider a family of distinct shift-invariant tessellations 
{Tj}j=o__j obtained by shifts of To. Their union covers W l J + 1-fold. We 
check refinability of the family, i.e. whether each member of the family can 
be expressed as a linear combination of the scaled-down copies of splines of 
the family. Example [2] makes this concrete for J = 2. 



, s := sin o- 











2s 




s/3 




'2s 2 ' 




3/2 


— s 




.-\/3/2_ 



Figure 3: A lozenge-shaped pair of triangles i is in the common support of 
three half-scaled translated copies of the grey hexagon h. 

Example 2 Denote by To a tessellation of the plane into unit-sized hexagons 
and by T\ and T2 its shifts by integer multiples of ^ni and 7^2 (see Fig. [3D . 
Let H E x(^o) be the indicator function of the unit hexagon h centered 
at the origin. Consider the three ^-scaled, translated copies of h shown in 
Fig. [3l The three copies intersect in a lozenge-shaped pair of triangles I. No 
other shifts of the ^-scaled hexagons in To, T\ or T2 overlap £. Therefore 
any linear combination of indicator functions in {Tj}j = o..j has a single value 
on I. Since i straddles the boundary of h, this constant linear combination 
would have to be simultaneously 1 and to replicate H. □ 

Example [5] suggests the following generalization of Proposition [5J A 
superposition T of a family of shift-invariant tessellations {Tj}j = Q j of R n 
is the tessellation obtained by partitioning W 1 by all cell facets in any of 
the tessellations Tj. T differs from {T,}j = o..j in that it contains fractions or 
pieces of the original cells. 

Proposition 4 Given a family {Tj}j = o..j of polyhedral shift- invariant tes- 
sellations ofM n , the space of indicator functions Uj=o jX(^i) * s refinable 
only if the superposition T of {Tj}j = o„j contains, for each facet f, the hy- 
perplane through f . 

Proof Assume that Uj=o jX(Tj) is refinable. Assume additionally that 
a scaled-down copy T 1 of T contains an unpartitioned piece c 1 E T 1 that 
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straddles a facet / of some cell c in one of the Tj. Analoguous to the 
proof of Proposition [TJ any linear combination of the indicator functions 
of the family {Tj}j = Q__j that replicates the step of the indicator function 
on c across / would have to be simultaneously and 1 on c . Therefore 
every facet of T must be a union of facets of T 1 . Analogous to the proof of 
Proposition [2j the claim follows by considering ever coarser-scaled copies of 
the Tj and hence of T. \\\ 
Example [3] illustrating Proposition [5] shows that Proposition 0] does not 
yield a sufficient constraint. In generalizing Proposition to overcomplete 
spline families, we restrict attention to families of tessellations that minimize 
facet overlap. 

Definition 1 (efficient family of tessellations) A family of tessellations 
{Tj}j=Q„j,Tj G W 1 is efficient if the intersection of more than two cell facets 
is of co-dimension greater than 1. 



Figure 4: (a) Two scaled-down shifts of a hexagon cover the top facet 
(edge) of the original hexagon exactly once from inside (gray region) and 
hence twice if we continue the tessellation by reflection across the facet. 
Adding the dashed scaled-down hexagon in (b) covers a part of the top 
(thick edge) twice from inside so that reflection yields an inefficient family 
of tessellations. 

Proposition 5 Let {2j"}j=o..J be an efficient family of tessellations ofW n 
by shifts of one polyhedral cell c G To. If two facets f a and /& of c meet with 
an obtuse angle a and no facet of c meets /& with an angle less than tt — a, 
and if d , the reflection of c across /&, is a cell of Tq then Uj = o jX(Tj) is 
not refinable. 

Proof Assume Uj=o jX(^j') * s refinable. Let e be the n — 2 dimensional 
intersection of two facets f a and /& of the cell c. By Proposition [H there 
exist c\,c\ G {Tj}j=o..j whose two facets / 16 and f\ h lie on /&, whose shared 
boundary e 1 is parallel to e (for n = 2, e 1 is a point) and that lie both to 
the same side of /& as c. Without loss of generality, f\ b is closer to e than 
is f\ h (cf. Fig. 2^,). When n = 2, we set p = e 1 and when n > 2, we pick a 
point p in the interior of e . 




(a) efficient 



(b) inefficient 
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Due to efficiency, of all cells in {Tj L }j = o..j that have one facet on fb and 
lie to the same side as c, exactly three pieces of the superposition T 1 of 
{Tj}j = o.j meet at p: the piece i± solely inside c\, the piece 12 solely inside 
C2, and in, the n-dimensional intersection of c\ and c?,. The intersection i n 
exists and is n-dimensional due to the obtuse angle a of c at e and hence 
of C2 at e ; and because c and hence c\ forms no angle less than ir — a with 
fb- By reflection across fb, there are three analogous pieces Oj outside c at 
p. Fig. [Sh, illustrates the situation in the 2-dimensional plane V through p 
and orthogonal to e . 




(a) 2 x 3 pieces (b) sliver inside in 




(c) sliver outside i n (d) overcomplete hex-splines 



Figure 5: Pieces (a) and slivers (red in b,c) in the plane V . (d) Concrete 
partition of Example [3] into pieces (sectors) of a neighborhood of point p on 
the partition of a facet /& of the support cell c of H. 

Denote by s the spline formed as a linear combination of indicator func- 
tions of the cells in \Tj }j=o..j that have one facet on To replicate 
the indicator function of c, s must have a unit step across fb and hence 
s(i\) — s(oi) = 1 and sfo) — s(o2) = 1 where s applied to a piece of T 1 
means evaluating s at some point in that piece, sufficiently close to p. Due 
to the overlap of the indicator functions on i n and on o n , at p 

s(i n ) - s(o n ) = s(h) + s(i 2 ) - s(oi) - s(o 2 ) = 2, (1) 

incompatible with the unit step across /& at p. 

Additional cells of {Tj} that overlap all six pieces, ii, o±, i n , o n , 12, 
02, surrounding p do not affect the above difference since their indicator 
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functions are constant. It remains to consider cells with facets crossing p and 
it suffices to consider the 2-dimensional plane V through p and orthogonal 
to e . In V, as shown in Fig. [5^,b,c, let the boundary (in, ^2) between i n 
and 12 form a smaller-or-equal angle with than the boundary (i n , Due 
to efficiency, within c, crossing facets lie either (b) strictly inside in or (c) 
strictly outside i n since the boundaries must not be covered a third time. 
In case (b) this strictness implies the existence of a sliver n, i.e. a piece of 
T 1 devoid of crossing facets, attached to (in ^2) and inside i n . Any linear 
combination of indicator functions corresponding to cells with crossing facets 
of type (b) therefore leaves unchanged the difference in value and hence the 
incompatibility (P), between and its reflection r D across (see Fig. [5b). 
In case (c), any linear combination of indicator functions that modifies the 
value in i n also modifies the difference in value between two slivers Gi and 
cr , where o\ is attached to (in,h) and lies inside i% and a Q is its reflection 
across ft,. Specifically, when the value of s(i n ) — s(o n ) is changed to satisfy 
the step condition, s(<Ji) — s(a a ) is changed away from the prescribed value 
1. Together this contradicts the assumption that Uj=o jX(Tj) is refinable. 

The following Example [3] illustrates Proposition [5l 
Example 3 Consider, as in Fig. [5]i, shifts 
Hi(x):=H(x-\-°]), H 2 (x):=H(x- 



7T 7T 

cos — , s := sin — 
3 3 



of the indicator function H{x) of a tessellation To. The three corresponding 
tessellations intersect only in single triangles and their superposition con- 
tains no piece that straddles the support hexagon of H. Yet this minimal 
familjH is, as expected, not refinable. Fig. [5b shows the regions in and 
o n where, for any spline s, s(i n ) — s(o n ) = 2 since s(ii) — s(oi) = 1 and 
s(i 2 ) ~ s(o 2 ) = 1. □ 

The proof of Corollary [2] showed that a pair of facets of the Voronoi 
cells of non-Cartesian crystallographic lattices form obtuse angles. By the 
symmetry of the cells, all facet angles are obtuse and hence satisfy the angle 
criteria of Proposition Together with the reflection symmetry of the 
lattices, Proposition [5] implies the following generalization of Corollary [2j 

Corollary 3 Efficient overcomplete families of splines obtained by convolu- 
tion of the Voronoi cell of a non- Cartesian crystallographic root lattice are 
not refinable. 



J = 2 yields the minimal number of scaled families that can satisfy the necessary 
constraints of Proposition [4] for hex-splines for three reasons. First, scaling should at least 
be binary. Second, the family that includes H/2 does not contribute to the step function 
because its cells straddle the boundary of c if they include a part of the boundary (see the 
dashed hexagon in Fig. [5ji). Third, at least two ^-scaled edges, namely of elements of Tl 
and Tl, are required to cover any edge of H. 
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5 Conclusion 



The paper identified several necessary criteria for tessellations to admit a re- 
finable space of (convolutions of) indicator functions. The criteria are chosen 
for their simplicity. For example, we showed that admissible shift-invariant 
tessellations must contain, for every facet, the whole plane through that 
facet. Already for hex-splines and hex-spline superpositions, an alternative 
algebraic proof of non-refinability is considerably more involved. 

Corollary [2] and [3] show that the increased isotropy of the Voronoi cells 
of non-Cartesian root lattices prevents refinability, even for overcomplete 
spaces obtained by efficient superposition of shifted lattices. Increased isotropy 
of the Voronoi cells is however the main reason for considering non-Cartesian 
lattices in the first place: these lattices have high packing densities that can 
improve sampling efficiency [PM62J. 

In conclusion, if we seek shift-invariant refinable classes of splines from 
convolving indicator functions of polyhedral cells, remarkably few options 
exist apart from tensor-product B-splines and box-splines. This does not 
imply that more general lattices fail to have associated refinable splines 
that represent their symmetry and translational structure. In the bivariate 
setting, odd orders of continuity on the hexagonal dual of the regular trian- 
gulation can be filled in by half-box splines [PB02] , And if fractal support 
is acceptable, [OS03, HR02J provide refinable functions with approximately 
hexagonal footprint. Combining families of symmetric box-splines, such 
as [KP11] yields refinable splines for any level of smoothness and crystal- 
lographic structure. It is just the particular approach of convolving non- 
Cartesian lattice Voronoi cells that fails to provide the important spline 
property of refinability. 

Acknowledgement Zhangjin Huang kindly worked out a first, alge- 
braic proof of non-refinability for Example [21 a scenario I posed to him. 
Andrew Vince and Carl de Boor helped me clarify the exposition in its early 
stages. 
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